The importance of local band effects for ferromagnetism in hole doped La2Cu04. 
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Band calculations for supercells of La(2_a;)Baa;Cu04 show that the rigid band model for doping is 
less adequate than what is commonly assumed. In particular, weak ferromagnetism (FM) can appear 
locally around clusters of high Ba concentration. The clustering is important at large dilution and 
averaged models for magnetism, such as the virtual crystal approximation, are unable to stabilize 
magnetic moments. These results give a support to the idea that weak FM can be the cause of the 
destruction of superconductivity at high hole doping. 
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The fascinating phase diagram of the cuprate high tem- 
perature superconductors still remains a great mystery in 
the field of condensed-matter physics |l| . Antiferromag- 
netic (AF) spin fluctuations are probably also relevant to 
the properties of the doped metallic state because they 
mix with phonons 0, S] and they provide mechanisms 
able to boost the electron-phonon interaction 4]. While 
many investigations have focused on underdoped phases 
fs, 6], little attention has been devoted to the overdoped 
regime, where the abrupt suppression of superconductiv- 
ity for a critical doping cannot be easily understood in 
term of the standard BCS theory. Only recently, Kopp 
et al. have suggested that ferromagnetic (FM) fluctu- 
ations compete with superconductivity in the overdoped 
regime. These authors have considered an itinerant pic- 
ture and have used it to justify the rich phenomenology 
of the overdoped cuprates. In particular, the termination 
of the superconductivity in the overdoped regime could 
mark a quantum critical point beyond which there should 
be a FM phase at zero temperature. A FM groundstate 
excludes a s or d wave superconductivity of the singlet 
variety The electron-phonon coupling, which may 
contribute to superconductivity, is also inhibited by FM 
spin fluctuations [1, [§| . 

Electronic structure calculations performed by Perry 
et al. ^Td\ using the spin unrestricted Becke-3-Lee- Yang- 
Parr (B3LYP) hybrid functional successfully produced a 
2.0 eV band gap and AF in (undoped) La2Cu04. How- 
ever, Paier et al. [l2| have raised serious concerns on 
the suitability of the B3LYP potential in metallic phases, 
and subsequent B3LYP calculations for minimal size su- 
percells of La2-a;Sra;Cu04, corresponding to x = 0.125, 
0.25, and 0.5, have predicted FM solutions. This 
in contrast to the non-magnetic or AF fluctuating states 
found in real, weakly doped systems. On the other hand, 
band calculations based on the local spin-density approx- 
imation (LSDA) yields the correct Fermi surface and a 
reasonable band dispersion and electronic structure 
studies on manganites have shown that the LSDA 
provides reliable potential for predicting both metallic 
and ferromagnetic phases. The difficulty with LSDA is 
rather to describe correctly the AF, insulating state of the 
undoped cuprates, but one should note that LSDA cal- 



culations for high-Tc cuprates with off-center lineariza- 
tion energies describe well both the undoped and the 
doped counterparts [l^. This fact indicates that two 
very different ground states are almost degenerate and 
that only small LSDA-corrections are needed to bring 
the two states into the correct order [T6|. 

In this letter, we confirm that the LSDA electronic 
structure of La2-2;Baa:Cu04 in the overdoped regime is 
consistent with weak ferromagnetism appearing locally 
around clusters of high Ba concentration. Our results 
from first principles validate the hypothesis by Kopp et 
al. [7\ that FM order and superconductivity are compet- 
ing in overdoped samples 

Our band calculations are performed with the Lin- 
ear Muflin-Tin Orbitial (LMTO) method [13, M, \M 
for supercells with two or 16 formula units (f.u.) of 
La2Cu04. The elementary 7-site unit cell is doubled 
along the x direction in the former case, which yields a 
(2,l,l)-extension of the elementary cell. The cell is then 
doubled along the 3 directions in the case of the large 
cell, leading to an (4,2,2)-extended cell with 112 atomic 
sites totally. The supercells in Ref. [ll| were adjusted 
to optimal size for each doping x, while here the use of 
a larger supercell of identical dimension for different x 
allow for some randomness in the local order. The basis 
set includes £ up through 3 for La- and Cu-sites, and up 
through 2 for 0-sites [20] . The number of fc-points is 84 
in the calculations for the small cell and it varies from 4 
to 125 in 1/4 of the Brillouin zone for the large cell. All 
sites in the cell are considered as non-equivalent through- 
out the self-consistent calculations. Other details of the 
methods can be found elsewhere [13, [S [H, [H HI] . 

Our result for paramagnetic La2Cu04 is in excellent 
agreement with the full potential calculation by Yu et 
al. [2^. As was mentioned above, doped La2Cu04 is 
metallic and LSDA calculations describe its band disper- 
sion and Fermi surface adequately, and AF can be pro- 
moted by use of off-set linearization energies in LMTO. 
The present calculations are based on a normal choice of 
linearization energies, and therefore we expect an under- 
estimation of the tendency towards magnetism, even in 
the doped cases. 

The magnetic moments for the three possible configu- 
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rations for a; = 1 in the supercells with two f.u.(14-site 
cells) are shown in Table IH Small but stable magnetic 
moments develop in the two configurations where the two 
Ba atoms in the cell belong to the same plane (i.e. are 
situated next to each other along x and y), or when they 
are shared between two planes. The two Cu in both 
configurations have almost equal magnetic moments and 
equal local density-of-states (DOS) at the Fermi energy 
Ep. Planar oxygen atoms have small magnetic moments, 
but the apical oxygen atoms between La/Ba and Cu has 
almost as large magnetic moment as Cu. Oxygen mag- 
netism is somewhat unusual, but there are cases reported 
in the literature like for a half-metallic Rb406 [13]. This 
participation of the apical O to the spontaneous magnetic 
polarization cannot be understood from the Stoner model 
because there is no clear correlation between local DOS at 
the apical O and the magnetic moment. Therefore, other 
mechanisms such as charge transfers and hybridization 
with neighbors are vital to explain the magnetization of 
the apical O atom [ll|. In our case, the charges on the 
apical O atoms are somewhat larger when these O sites 
are between the Cu atom and the Ba atom rather than 
between Cu atom and La atom, while the magnetic mo- 
ments are largest on the O atom close to the La-plane 
when the La/Ba atoms are fully occupying a plane. This 
trend is reversed when there is alternate La/Ba occupa- 
tion within the planes. No magnetic moment is developed 
when the two Ba atoms are located on top of each other, 
i.e. along z. This last result might seem surprising, since 
a finite magnetic moment is present next to a single Ba 
atom in the two other cases. 

Therefore, it seems that inhomogeneous formation of 
doped holes in the vicinity of the impurity and the clus- 
tering of impurity atoms within the same plane are im- 
portant ingredients for the formation of magnetic mo- 
ments on Cu atoms and apical O sites. This conclusion 
is corroborated by the fact that no magnetism is found 
when the virtual crystal approximation (VCA) is applied 
for X = 1. The missing charge on Ba sites relative to La 
sites is only partly taken from the Cu band, so Ep does 
not quite reach the high DOS region, thus the system is 
below the Stoner limit for FM. Eventually, when the Ba 
concentration becomes larger (x > 1.5), FM appears in 
VCA calculations. Random oxygen vacancies, described 
by CPA calculations, will increase the DOS at Ep [25j . 
and might help Stoner magnetism. 

Next, by using larger supercell calculations with Ba/La 
substitutions, we try to answer the question if local spon- 
taneous magnetization is possible for hole doping corre- 
sponding to the most "over-doped" (actually meaning 
low doping for supercell calculations) high-Tc materials, 
and thus if FM can be responsible for the disappearance 
of superconductivity. For the large supercells with 16 f.u. 
(112-site cells), we present results for configurations with 
4, 10 and 16 Ba sites distributed randomly or with a ten- 
dency of clustering on the 32 available La sites. We focus 
on the clustered configurations since they show enhanced 
spontaneous magnetization. The calculations with 4 Ba 



sites correspond to the interesting doping concentration 
X = 0.25. 

We first consider the paramagnetic state for 4 and 10 
clustered Ba atoms, where the DOS are calculated within 
a larger mesh of 125 k points. While the main DOS fea- 
tures with random distribution of Ba sites are close to 
a rigid band picture given by our VCA calculations and 
by the CPA results for a; = 0.14 |25i] and for x = 0.30 
[20 |. the DOS for clustered configurations can be differ- 
ent. The results for the clustered cases are compared in 
Fig.[l]together with DOS for the pristine phase La2Cu04 
with a rigid band shift of Ep. As shown in Fig. (TJ even 
a low concentration of Ba atoms is sufhcient to induce 
band broadening and shifts of the DOS peaks towards 
Ep, whereas the rigid band model only gives a small 
DOS value at Ep. Broadening due to dilute substitu- 
tion of La for Ba is not the whole story since the local 
DOS on individual Cu sites can also change dramatically 
during the self-consistent procedure of the spin-polarized 
calculations. Curiously, the spin-polarized calculations 
show that the largest local magnetic moments and their 
corresponding exchange splitting of the order of 1 mRy 
do not always occur on sites with large local paramag- 
netic DOS. Therefore, while the Stoner criterion seems 
to provide an overall indicator for possible magnetiza- 
tion, the paramagnetic DOS may be different from the 
actual shapes of majority and minority DOS in the FM 
phase. These strong non-linear effects are also reflected 
by a slow convergence of the magnetic moment during the 
self-consistent procedure. Therefore, the results for the 
large supercell confirm that the charge transfer from the 
dopants to the CuOg octahedra adjacent to the Ba impu- 
rities and the resulting spontaneous magnetization can- 
not be modeled by using rigid band and Stoner schemes. 

For the spin-polarized calculations, we also have stud- 
ied the behavior of the magnetic moment as function of 
the number of fc-points. The largest magnetic moments 
occur for simulations with only 4 fc-points. These calcu- 
lations are analogous to models for FM nano-domainsj? | , 
where the pinning of Ep at DOS peaks is observed [27 j. 
However, in order to simulate bulk properties, we have 
considered computations with 12, 27 and 125 /c-points. 
The total magnetic moments as function of fc-points and 
different Ba configurations are displayed in Tabic [TTl The 
magnetic moments are localized on the Cu and apical O 
sites which are close to Ba sites, as in the case for the 



TABLE I: The calculated magnetic moments (/is per cell) for 
La2Ba2Cu208 for the 3 configurations of Ba positions. The 
occupations of Ba/La are denoted as "pl-l" for plane-1, "col- 
2" for the second column, etc... 
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FIG. 1: (Color online) Total Paramagnetic DOS for the 
La28Ba4Cui6064 and the La22BaioCui6064 clusters com- 
pared to a rigid band model for La2Cu04, in which Ep is 
shifted to correspond to x — 0.62. 



small cell. This can be seen on Fig.[2]and Fig.[3]of the su- 
percell where the lengths of the arrows are proportional 
to the local magnetic moment on each site |28j . 

As shown in Fig.[2l both the 16 Ba-site supercells, ran- 
dom and cluster, yield magnetic solutions. Therefore, a 
random distribution of Ba sites does not destroy mag- 
netic solutions for very high Ba doping. Fig. [5] confirms 
that the apical O sites close to Ba impurities have the 
highest magnetic moments. 

Fig. [3] illustrates the appearance of spontaneous mag- 
netic moments for lower Ba concentrations. When de- 
creasing Ba concentration to 10 clustered atoms, the api- 
cal O sites do not develop large magnetic moments as for 
the configurations with 16 Ba atoms. The effect of clus- 
tering is also seen in the case with four Ba atoms, i.e. 
for doping corresponding to the interesting range of dop- 
ing for "overdoped" superconducting samples. With four 
clustered Ba atoms in the 112-site cell (3 next each other 
and the 4th one opposite in the next layer) it is possible 



TABLE II: The calculated magnetic moments (hb per cell) for 
La(2-x)Baa;Cu208 as function of fc-points for 5 configurations 
of Ba substitution with 16, 10 and 4 Ba on the 32 possible La 
positions. The averaged separation between energy levels is 
lower than the exchange splitting for the most dense fc-point 
mesh, which indicates a sufficient fc-convergence. 



FIG. 2: (Color online) Magnetic moments (arrows pro- 
portional to the size of the moments) on the atoms of the 
Lai6Bai6Cui6064 (a) random and (b) cluster supercells (only 
non-equivalent atoms are shown). The (gray) small spheres 
are Ba atoms, the small (blue) spheres are La atoms, the large 
(green) spheres are Cu atoms and the smallest (red) spheres 
are O atoms. The limits of an octahedron are outlined by yel- 
low in (b). The calculations have been performed with 27-k 
points. 
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FIG. 3: (Color online) Same as Fig. [2] for the (a) 
La22BaioCui6064 and (b) La28Ba4Cui6 064 supercells. The 
arrows have been rescaled by a factor of 2 compared to Fig. [2] 
for clarity. 
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to follow a contamination of the induced magnetic mo- 
ment towards the Cu atoms in the second layer as shown 
in Fig. [Sljb) . About 60 % of the total magnetic moment 
is on this Cu plane, which has less electrons than the first 
one. Interestingly, this self-consistent result has been ob- 



tained from a starting configuration where the magnetic 
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moments were put near the Ba in the first layer of Cu. 
The magnetic moments on the O atoms (apical and pla- 
nar) are generally very small at this low doping. This is in 
contrast to the result at larger doping, where apical sites 
could act as " links" between magnetic layers and thereby 
acquire magnetic moments. When the four Ba atoms are 
distributed randomly (i.e. not next to each other) the 
magnetic moment goes to zero (cf. Table HI)) . This un- 
derlines the necessity for having a rather strong, local 
perturbation (a few Ba atoms within the same plane) in 
order to polarize neighboring Cu atoms. 

Our results agree with the experimental findings that 
magnetization starts near x = 0.22 and corroborate the 
interpretation of neutron scattering data in terms of in- 
homogeneities such as Ba clusters and FM nano-domains 
[29i] . The observed large increase of the magnetic moment 
with doping [2^ may indicate a substantial clustering. 
In our calculations, the highest magnetic moment for al- 
ternate La/Ba layers at x—1 (Table [l| corresponds to 
m — I.Ghb per large cell instead of m = 0.43/LtB for the 
random configuration in Table [U [sO] ■ Furthermore, the 
case with x = 0.25 and 4 fc-points, which can simulate 
an isolated FM nano-domain as mentioned above, gives 
a magnetic moment of to = 0.82/iB which is close to the 
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the doping is high, and even at a; = 0.25 if there is some 
clustering of the Ba atoms. In the latter case, only some 
regions near the Ba clusters are weakly polarized. The 
mechanism of FM from Ba- (or Sr-) clusters could be 
responsible for the suppression of superconductivity in 
"overdoped" cuprates as suggested by Kopp et aZ.Q- 
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